Majorana zero modes in a quantum Ising chain with longer-ranged interactions 
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A one-dimensional Ising model in a transverse field can be mapped onto a system of spinless fermions with 
p-wave superconductivity. In the weak-coupling BCS regime, it exhibits a zero energy Majorana mode at each 
end of the chain. Here, we consider a variation of the model, which represents a superconductor with longer 
ranged kinetic energy and pairing amplitudes, as is likely to occur in more realistic systems. It possesses a richer 
zero temperature phase diagram and has several quantum phase transitions. From an exact solution of the model 
these phases can be classified according to the number of Majorana zero modes of an open chain: 0, 1, or 2 at 
each end. The model posseses a multicritical point where phases with 0, 1, and 2 Majorana end modes meet. 
The number of Majorana modes at each end of the chain is identical to the topological winding number of the 
Anderson's pseudospin vector that describes the BCS Hamiltonian. The topological classification of the phases 
requires a unitary time-reversal symmetry to be present. When this symmetry is broken, only the number of 
Majorana end modes modulo 2 can be used to distinguish two phases. In one of the regimes, the wave functions 
of the two phase shifted Majorana zero modes decays exponentially in space but in an oscillatory manner. The 
wavelength of oscillation is identical to the asymptotic connected spin-spin correlation of the XF-model in a 
transverse field to which our model is dual. 



I. INTRODUCTION 

There has been much recent interest in Majorana zero 
modes.EEII Their relevance to topologically protected quan- 
tum computation is intensely studied. Kitaev '"^ suggested an 
elegant model of a one-dimensional p-wave superconducting 
wire, which supports Majorana zero modes at the ends of the 
chain. 

Kitaev's model is the fermionized version of the famil- 
iar one-dimensional transverse field Ising model (TFIM),'^ 
which is one of the simplest models of quantum critical- 
ity. In the fermionic representation, the well-known quantum 
phase transition in the model can be understood as a transition 
from the weak pairing BCS regime to the strong pairing BEC 
regime. The weak-pairing phase is topologically non-trivial 
and in this phase the chain with open boundaries possesses a 
Majorana fermion zero enegy mode localized at each end. It 
is equivalent to the ferromagnetic phase of the transverse field 
Ising chain. The strong-pairing phase is topologically trivial 
and does not have any normalizable Majorana fermion zero 
enegy modes at the ends. It corresponds to the quantum disor- 
dered phase of the transverse field Ising chain. Recently, there 
have been attempts to realize Kitaev's model in one dimen- 
sional wire networks.'^ In a realistic quantum wire, however, 
the range of the hybridization of the electron wave function, as 
well as that of Cooper paiiing will be of finite range and the 
effect of such longer ranged interactions must be addressed. 
The goal of this paper is to study the effect of such longer 
ranged interactions. We do so within the context of another 
exactly solvable model, and find a rich phase diagram that re- 
sults from such longer ranged correlations. 

These longer ranged interactions were considered in a pre- 
viously introduced generalization of the Ising model in a 
transverse field by extending it to contain a three spin interac- 
tion term,n3 which is also exactly solved by a Jordan- Wigner 
transformation.n^This generalized model arises as a first step 



of a real space renormalization group transformatiorP^and has 
a richer phase diagram. The purpose of that study was to un- 
derstand how irrelevant operators can drive a system along a 
critical line between two different zero temperature quantum 
critical points. The flow of this crossover as the higher energy 
states are integrated out conform to the Zamolodchikov's c- 
theorem.^" It is an explicit example of an exactly solved case 
where a system that points to a given fixed point at higher 
energies can asymptotically flow to a different fixed point at 
lower energy scales. This flow was explicitly traced in terms 
of a flow from higher temperature to lower temperature. The 
lesson learnt there was that at higher temperatures a system 
may be pointing to a different fixed point compared to its true 
fate at zero temperature. 

Here we reexamine the Ising model with a transverse field, 
with the added three spin interaction from the perspective 
of Majorana zero modes. We find that the phase diagram 
can be classified according to the number of Majorana zero 
modes. The fermionized version of this model corresponds to 
a p-wave superconductor in which the electrons have longer 
ranged hoppings and longer ranged harmonics of the p-wave 
gap function, enabling us to address the effect of such longer 
ranged interactions on the zero temperature phase diagram of 
the quantum Ising chain. We find several topological phase 
transitions in our model, and the phases can be classified by a 
topological invariant of the Anderson pseudospin vectoi^ of 
the mean-field description of the superconducting state. This 
topological invariant is an integer, Z, and also specifies the 
number of normalizable Majorana fermion zero energy modes 
that are localized at each end of a chain with open boundary 
conditions. 

We assess the conditions under which the topological or- 
der of the zero temperature phase diagram remains intact, and 
find that all phases are protected by a unitary version of time- 
reversal symmetry (appropriate for spinless fermions): so long 
as this time reversal symmetry is preserved, the phases de- 
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scribed in our work remain stable. In particular, we find that 
even when there are 2 Majorana fermion zero modes localized 
at each end, separated by a lattice spacing with wave functions 
orthogonal to each other. Once we allow breaking of time 
reversal invariance the topological invariant collapses to Z2, 
which implies at most one Majorana zero mode at each end of 
an open chain. The results of the ID superconductor also of- 
fers insights into TFIM. We find that the phase dominated by 
the three spin interaction has ground state degeneracy from 
analyzing the Majorana zero modes. However, we also find 
that there exists a class of local spin interaction that can re- 
move this ground state degeneracy. Such impurities are time 
reversal breaking and it is not clear how they realize in generic 
circumstances. 

The dual^^ (exchanging the site spins by the bond spins) 
of the three-spin model that we study is amusingly the one- 
dimensional quantum XF-model in a transverse magnetic 
field^^ in an enlarged parameter space than studied previously. 
From a complex calculation of the asymptotic form of the 
connected z-component of the instantaneous spin-spin corre- 
lation function it was discovered that there is an oscillatory 
region within the ferromagnetic phase. We find that this phe- 
nomenon of oscillation is intimately related to the oscillation 
of the Majorana zero modes, most remarkably the oscillation 
wavelengths are identical, as is the exponential decay in the 
vicinity of the quantum critical lines. 

The plan of the paper is as follows: in Sec. II we set the 
stage by recapitulating the phase diagram of the model to ori- 
ent the reader In Sec. Ill Majorana zero modes and their 
properties are obtained from the solution of the Bogoliubov- 
de Gennes (BdG) equation, while in Sec. IV we discuss the 
efficacy of the Majorana representation by obtaining the so- 
lution of a three-term recursion relation instead of the full 
numerical solution of the BdG Hamiltonian. In Sec. V we 
discuss the topological aspects and Sec. VI is the concluding 
section. There are three appendices giving some details. 



II. THE HAMILTONIAN AND THE PHASE DIAGRAM 
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H = -5^1(1 - 2c|c,) - Ai ^ (4c,+i + 44+1 + h.c.) 

i=l 1=1 
- A2 ^ (4-lCj+l + Ci+lQ_i + /i.e.). 



1=2 



(4) 



In contrast to the spin model, the spinless fermion Hamil- 
tonian is actually a one-dimensional mean field model for a 
triplet superconductor, where there are both nearest and next 
nearest neighbor hopping, as well as condensates. The near- 
est neighbor hopping amplitude Ai is also the amplitude of 
the nearest neighbor superconducting gap, and the next near- 
est neighbor hopping amplitude is equal to the next nearest 
neighbor superconducting gap; in general Ai 7^ A2. In terms 
of Jordan-Wigner fermions one can envision finding an actual 
one-dimensional system with such an extended Hamiltonian. 
The solution of the corresponding spin Hamiltonian through 
Jordan-Wigner transformation is, however, exact and includes 
all possible fluctuation effects and is not a mean field solution. 

Imposing periodic boundary condition, the Hamiltonian 
can be immediately diagonalized by a Bogoliubov transfor- 
mation: 
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(5) 



The anticommuting fermion operators jy^'s are suitable linear 
combinations in the momentum space of the original Jordan- 
Wigner fermions. The spectra of excitations are (lattice spac- 
ing will be set to unity throughout the paper) 

efe = ±2^ 1 + Af + A^ + 2Ai(l - A2) cos k - 2X2 cos 2k 

(6) 

unless, otherwise stated, we shall set g = 1. Quantum phase 
transitions of this model are given by the nonanalyticities of 
the ground state energy: 



^0 - ^Ty^^k- 



(7) 



The three spin extension of the TFIM, which was previ- 
ously studied,''^ has the Hamiltonian 

H^- Y.^gcJ- + A2af af_iaf+i + \i(Jtal_^) (1) 



The cr's are the standard Pauli matrices. In this section we in- 
troduce the Hamiltonian and its phase diagram from a conven- 
tional Jordan-Wigner analysis. The Hamiltonian after Jordan- 
Wigner transformation 
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11(1 -2c]c,)(c,+ 4) 



(2) 
(3) 



These nonanalyticites are also defined by the critical lines 
where the gaps collapse; see Fig.[T] 

For the Ising model in a transverse field without three spin 
interaction, the gaps collapse at the Brillouin zone boundaries, 
k — ±7r at the self-dual point Ai = 1 and A2 — 0. When the 
three spin interaction is added, the gaps can collapse at fc = 
as well as fc = arccos(Ai/2) for A2 = —1 and < Ai < 
2. But at the free fermion point b, there are no zero energy 
excitations except at fc = ±7r/2. When we move to the point 
c the spectrum evolves increasing the weight at fc = and the 
locations of the nodes are incommensurate with the lattice. 
The incommensuration shifts as a function of A2 . The point d 
is a multicritical point and the spectra vanishes quadratically 
at ±7r. As we shall discuss below, the spectra are no longer 
relativistic at this point as a result of the confluence of two 
Dirac points, corresponding to a dynamical exponent z = 2. 
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FIG. 1. (Color online)The region A2 > 1 + Ai is ordered in the 
original spin representation and the boundary of it is a critical line 
where the gap at fc = collapses. The region A2 < 1 — Ai is disor- 
dered as well and the boundary corresponds a critical line where the 
gap at fc = TT collapses. The point Ai = and A2 = 1 is a special 
multicrtical point with an emergent (7(1) symmetry, most transpar- 
ently seen in the dual representation (see below). In the same dual 
representation, the region enclosed by Ai = — 4A2 is an oscillatory 
ferromagnetically ordered phase separating from an ordered phase 
for A2 < 0, as determined by the spatial decay of the instantaneous 
spin-spin correlation function. Note that duality exchanges ordered 
and disordered phases. Here n — 0,1,2 correspond to regions with 
n-Majorana zero modes at each end of an open chain. 



III. MAJORANA ZERO MODES 



open boundary condition. The phase diagram itself can be de- 
duced from the number of zero modes of the BdG equation. 
In the next section we shall see that in the Majorana repre- 
sentation, the zero modes can be obtained from a very simple 
recursion relation. With reference to Fig. [T] we note that there 
are regions of n = 0, n = 1, and n = 2 zero modes, and 
the lines separating them are quantum critical lines, except 
for the line separating n ~ and n = 2, which is a topolog- 
ical transition. The thin line A^ = — 4A2 corresponds to zero 
entanglement entropy. In this respect, it is remarkable that 
this thin line osculates the quantum critical line. 

For A2 > 0, and A2 > 1 + Ai, one of the zero modes decays 
exponentially in the bulk, and the decay length diverges as 
the quantum critical line is approached. The amplitude of the 
second zero mode also decays exponentially but it oscillates 
as e^™ regardless of Ai as it approaches the quantum critical 
line, at which point it loses the compactness of its support, 
signifying the loss of this zero mode. On the side A2 < 1 + Ai, 
one zero mode is recovered and it decays exponentially in the 
bulk, as in the region A2 > 1 + Ai, as shown in Fig.|2] 
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In this section we explore the zero modes by the 
Bogoliubov-de Gennes equations with open boundary condi- 
tion. 



A. Unbroken time reversal invariance 

The equations, assuming open boundary condition, are 
given by 



h A 
-A -h 



(8) 



where the submatrices are (unless otherwise stated, we will 
set 3 = 1) 

hij = Ai((5j^i+i + Sjs^i) + A2((5j,i+2 + - 2%(9) 

\j = -Ai((5j,i+i - - A2((5jVi+2 - (10) 

Here — (u„(l), u„(2), . . . u{N)) and from time reversal 
symmetry of the Hamiltonian m„ = m* and Vn — v^- The 
eigenvalue is labeled by n and the arguments of u and v are 
lattice indices. 

From the diagonaliztion of the BdG Hamiltonian, we can 
easily see that the Majorana zero-modes can occur only for 
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FIG. 2. (Color online)The two Majorana zero modes for Ai — 0.05 
and A2 — 1.5. The row number 1 shows Moi(i) and voi{i) cor- 
responding to the the first Majorana zero mode. The second row 
corresponds to the second Majorana zero mode that is orthogonal to 
the first. The third row corresponds to the probability distributions 
Poi{i) and Po2(*) for the two respective Majorana modes. The nu- 
merical diagonalization was carried out for a lattice of TV = 100 
sites with open boundary condition. For lattices larger than 150, one 
quickly looses numerical control because of 10^'' difference in the 
order of magnitudes of the largest eigenvalue and the zero mode. 

The situation is richer for A2 < 0. First of all there are 
no zero modes until A2 < — 1 for A2 < 1 — Ai. In this 
region, there are two zero modes, both of which are oscilla- 
tory with exponentially decaying envelope. But this time the 
wavelengths of the modes depend on the parameters (Ai, A2). 
Note that they are phase shifted by a lattice site; see Fig. |3] 
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When we cross the quantum critical line, A2 = 1 — Ai, a non- 
oscillatory and exponentially decaying zero mode is observed. 
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A2 will be purely imaginary, while the Ai term will remain 
untouched. This is easier to compare with the previous phase 
diagram because the absolute magnitude of the next-nearest 
neighbor condensate remains the same. We have verified that 
our conclusions hold for arbitrary 9 as well. 





FIG. 4. (Color onIine)Splitting of n — 2 Majorana zero modes for 
the complex Hermitian BdG equation. Left: A2 — 2.5. Right: A2 = 
— 2.0. The magnitude of eo is the smallest eigenvalue. The slight 
rounding in the proximity of the quantum critical point is due to the 
finite size of the lattice: A'^ — 200. The quantum critical point for the 
infinite system is at Ai = 1.5 for the left panel and Ai — 3.0 for the 
right panel. The n = 1 Majorana zero mode survives intact. Note 
that at Ai =0 the chain splits into two independent chains and hence 
there is a zero mode irrespective of the fact that A2 pairing amplitude 
is purely imaginary. 



FIG. 3. (Color online)The two Majorana zero modes for s A2 = 
— 1.2 and Ai = 1. The row number 1 depicts uo2(j) and vo2{i) 
corresponding to the the first Majorana zero mode. The second row 
shows the second Majorana zero mode that is orthogonal to the first. 
The third row corresponds to the probability distributions Poi (i) and 
P02 (i) for the two respective Majorana modes. 



B. Broken time reversal invariance 

We can ask what happens if we add a relative phase between 
the two order parameters in the BdG Hamiltonian, while keep- 
ing the single particle Hamiltonian intact. Then, 



h A 
At -h 



where the submatrices are 



(11) 



- 26,, (12) 
Ay = e*''A2(5j,j+2 - Si+2,j) + Ai(*i_ij - (5i+ij)(13) 

The solution of this modified BdG Hamiltonian shows that 
the regions of the phase diagram which contain n = 1 Majo- 
rana zero mode remain robust while those containing n ~ 2 
Majorana zero modes are in general destroyed, meaning that 
they are split; two examples are shown below in Fig. Q for 
6 = Tr/2. For an arbitrary value of 9, the real and the imagi- 
nary parts of Aij will receive contributions from both Ai and 
A2 making it difficult to directly compare with the previous 
phase diagram in Fig.[T] For 9 ~ 7r/2, the term containing 



IV. MAJORANA HAMILTONIAN 

To unveil the hidden topology behind the TFIM with three 
spin interaction, we introduce Majorana fermion operators: 



(14) 
(15) 



satisfying anti-commutation relations, so that the Hamiltonian 
in Eq. (j4]) becomes 



H 



N 



iV-1 



i=l 



i=2 



(16) 

The three-spin interaction corresponds to a next-nearest 
neighbor coupling in addition to the Majorana fermion ana- 
log of polyacetylene in the transverse ising model. The differ- 
ence here is that onsite-potentials do not occur with Majorana 
fermions: terms like Va^ or Vb^ only add overall constants, 
since = 6^ = 1 We write down the Hamiltonian in the 
basis ■i/;^ = (oi, 61, 02, 62, i3: ^3---) and see that 
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(17) 



The zero modes are given by the recursion relation of the 
amplitudes. 



A, + \iAi+i + A2A,+2 = 0, 



(18) 
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where the eigenvector is chosen to be of the form 
(^1, 0, ^2, 0, • • • )^. We note that when A2 < 0, the above re- 
cursion relation looks like the equation of motion of a damped 
harmonic oscillator with the time variable discretized. The 
two linearly independent solutions can be expressed as: 



and 



where q± satisfy 



1 = \iq + \-2,(^ 



<1± 



2A2 



(19) 

(20) 
(21) 



and Ci, C2 are constants. However, when we restrict our- 
selves to real and normalizable solutions, we may have only 
one, two, or zero solutions. Note that + 4A2 > is suffi- 
cient for giving us real solution regardless of Ci to C2 ratio. 
Therefore, we can obtain the phase boundary for Af +4A2 > 
just by examining whether \q± \ is larger than 1 or not. For the 
case Ai, A2 > 0, we can recover the results found by Kopp 
and Chakravarty. For 1 — Ai < A2 < 1 + Ai, we have 
< (7+ < 1 < — g_, and therefore a single Majorana zero 
mode at each end of the chain. For A2 > 1 + Ai, we have 
< — < 1 and < (7+ < 1 and thus there are two Majo- 
rana zero modes at each end of the chain. Lastly, we have no 
Majorana zero mode for A2 < 1 — Ai, as (7+ and —q_ are both 
greater than 1 . We can extend this analysis to Af + 4A2 > 
and A2 < 0, Ai > 0. For A2 < 1 - Ai and A2 > -1, 
we have \q±\ > 1 and thus no Majorana zero mode, whereas 
A2 > 1 — Ai gives us < < 1 < — g_ and thus there is 
one Majorana zero mode at each end of the chain. 

By contrast, when we break the time reversal symmetry by 
having phase difference between the two pairing terms as in 
Eq. (13 I, we find that we can have only zero or one normal- 



izable Majorana zero mode at the end of the chain, as shown 
above from the explicit solution of the BdG equation. Interest- 
ingly, we find that we obtain the same result if we break time 
reversal symmetry by adding an impurity term of the form 
(see Appendix [B) 



j-\-7n 



(22) 



to the original Hamiltonian Eq. Q. Since the translational 
invariance is broken in this case, the number of zero Majo- 
rana modes provide a convenient way to distinguish different 
phases. 



A. Oscillatory Majorana zero modes with varying wavelength 



We find from the recursion relation, Eq. ( frsj l, that there are 
no Majorana zero modes for — 1 < A2 < and A2 < 1 — Ai. 
However, there are two oscillatory zero modes for A2 < — 1 
and A^ + 4A2 < 0, with amplitudes at a lattice site j given by 
the two solutions of the recursion relation: 



A, 



(-A2) 



-3/2 



cos J 6 



(23) 



A, = (-A2)--'/2sinj0, 



(24) 



where 9 = arcsin(Ai/-\/— 4A2). The amplitude could be 
rewritten as 



(-A2) 



(25) 



where ^ = 2a/ ln(— A2). Note that close to the quantum crit- 
ical line IA2I ^ |Ai — 1|. We have reintroduced the lattice 
spacing a here. An example of oscillatory Majorana modes 
are shown in Fig. ([3]). 

When A2 becomes negative, an oscillatory phase, as deter- 
mined from the spin-spin correlation function, was obtained 
from a dual transformation that exchanges sites and bonds of 
the lattice. Then the Hamiltonian in Eq. Q can be cast in 
the standard notation of the quantum XF-model by factoring 
out an overall scale. Thus, with /x's as bond-centered Pauli 
matrices. 



H 



1 + r 



1 + r 



fii{n)^i{n + 1) 



1 - r 



^J■2in)^l2{n + 1) + h^j.3{n)) 



(26) 



the two parametrizations are related to each other by 
2h 



Ai 



1 + r 



A. 



r- 1 
1 + r ■ 



(27) 



The critical line in XY-model, separating the quantum dis- 
ordered phase from the ferromagnetic phase, is ft, = 1, which 
corresponds to Ai + A2 = 1, separating the ordered phase from 
the disordered phase. The model was previously studied only 
in the range < r < 1 and h > 0. Since the ordered and the 
disordered phases are exchanged under duality, the disordered 
phase of the three-spin model is Ai + A2 < 1. 

A complex calculatiorPlof the instantaneous spin-spin cor- 
relation function showed that within the ferromagnetic phase, 
there is an oscillatory phase in the which the connected corre- 
lation function G(a::) = {^3{x)fi{0)) — {firj{x)) {fi3{Q)) in the 
limit a; 00, is 



G{r) 



1 <=-^/« 



disordered, 
ordered, 

oscillatory ordered. 



(28) 



Here cosK — Ai/V~4A2 and f is the spin-spin correlation 
length. The oscillatory phase in the XF-model is bounded by 
r"^ + h? < 1, which corresponds to A2 < — A^/4 in the three- 
spin model. Note that the oscillation wavelength is identical 
to the wavelength of the Majorana fermions. Even the corre- 
lation length close to criticality is the scale of the exponential 
decay of the Majorana fermions. This must imply that in the 
spectral decomposition, the Majorana zero modes asymptot- 
ically dominate, although we have not yet found a rigorous 
proof of it. 

Since Majorana modes of zero energy are degenerate eigen- 
states, a change of the number of Majorana zero modes can 
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only occur when the energy gap collapses, i.e. at a quan- 
tum phase transition. Reexamining the behavior of Majorana 
zero mode in the full parameter space, there are three dividing 
lines based on their number: A2 = 1 + Ai, A2 = 1 — Ai and 
A2 — —1; see Fig[T] These lines are identified with the critical 
lines signifying phase transition, as can be seen in the energy 
spectrum in the previous result. So in this case, the number 
of Majorana zero modes serve as an "order parameter" for the 
quantum phase transition. 

We can then distinguish phases and locate quantum phase 
transitions by simply transforming the Hamiltonian in terms 
of Majorana operators and finding the number of allowed Ma- 
jorana zero-modes at each end of the chain, which is a topo- 
logically protected quantity. This provides us a profoundly 
simple way to study quantum phase transitions. 

B. Unbroken unitary time reversal symmetry 

Some insight into the phase diagram can be obtained from 
the perspective of the weak to strong pairing topological phase 
transitions in this model. Starting from the spinless Fermion 
Hamiltonian in Eq. (j4]), we get upon Fourier transformation 

H = (2 — 2Ai cos fc — 2A2 cos 2A:) cJ,Cfc 

/ \ ^^^^ 

+ f lAi sin fccj^c^ J, + zA2 sin 2fcc[,c^ j, + /i.e. j , 

fe 

which describes a superconductor with a pairing potential that 
consists of a nearest neighbor and a second nearest neighbor 
p-wave pairing. The BdG Hamiltonian which governs the dy- 
namics of the BCS quasiparticles at each momentum fc has the 
form 

= ( -iA(S r-";:. ) <»> 

where = — 2AiCosfc — 2A2COs2A:, A(A:) = Aisinfc + 
A2 sin2fc, and /i = —2. In this representation of the model, 
the various phase boundaries described in previous sections 
correspond to Lifshitz transitions, across which the number 
of Fermi points change. However, as we shall see, not all 
Lifshitz transitions are topological phase transitions. To de- 
termine precisely whether a Lifshitz transition is a topological 
phase transition, we must define an integer-valued topological 
invariant that changes only across a topological phase transi- 
tion. It is convenient to define the invariant using the Ander- 
son pseudospin vectoi^^ 

d(fc) = A(fc)y + (£,.- (31) 

In terms of this vector Hsdcik) = d{k) ■ t, where f are 
Pauli matrices which act in the Nambu (i.e. particle-hole) ba- 
sis of HsdG- It is important to highlight that the pseudospin 
is defined only in the yz plane in this problem. This is a 
consequence of time-reversal symmetry (applied to spinless 
Fermions, time-reversal is simply the operation of complex- 
conjugation). Time-reversal symmetry ensures that the rel- 
ative phase between the nearest neighbor and second neigh- 
bor pairing amplitudes must be real (in the context of the spin 



model, the relative phase is identically zero). However, when 
time-reversal symmetry is broken, the relative phase between 
these can be an arbitrary complex number. If this happens, 
the Anderson pseudospin vector will have three components, 
and the analysis below is invalidated. In this section, we shall 
restrict our attention to the case where the relative phase is 
zero. The topological invariant that characterizes the phase 
transitions will be defined in terms of the unit vector 

dik) 

d{k) = — = cos Okij + sin^fcz. (32) 

\d{k)\ 

Here, the momentum states with periodic boundary conditions 
form a ring T^, and the unit vector d{k), lives on a unit circle 
in the yz— plane. Therefore, the angle 6{k) is a mapping 
9{k) : —>■ and the topological invariant we seek is 
simply the fundamental group of this mapping, which is just 
the integer winding number 

W = j^ (33) 

where the integral is done around the one dimensional Bril- 
louin zone. This quantity characterizes the number of times 
the vector d{k) rotates in the j/z-plane around the one- 
dimensional Brillouin zone. It can only be an integer and 
therefore cannot vary with smooth deformations of the Hamil- 
tonian, so long as the quasiparticle gap remains finite. The 
winding number changes discontinuously only when the en- 
ergy gap vanishes, i.e. at a topological phase transition. More- 
over, the change in the number of normalizable Majorana 
modes at each end of the chain across a transition is given 
by the change in the winding number We now apply this 
framework to characterize several critical points in the A1A2 
plane. In Fig.|5]we show the results in the vicinity of the criti- 
cal point at Ai = 1 and A2 = 0. At this critical point, the chain 
consists only of nearest neighbor hopping and pairing. There- 
fore, the kinetic energy has its minimum at fc = 0, and the gap 
function also vanishes at this point. For Ai = 1^ and A2 = 0, 
the chemical potential occurs below the band bottom. In this 
limit, the winding number of d{k) is zero, since the configura- 
tion is topologically equivalent to one where 9{k) — 7r/2 for 
all fc, and this state clearly has W — Q. On the other side of 
the transition, Ai = 1+ and A2 = 0, the chemical potential 
crosses the band bottom, and the winding number changes to 
W — 1. At the critical point itself, the Anderson pseudospin 
unit-vector is not defined at the point fc = 0, where the gap 
closure occurs. The change in W is identical to the change 
in the number of normalizable Majorana fermion zero modes 
across this transition. 

Next, we use similar reasoning to study the transition at the 
critical point A2 = 1 and Ai = across which the change 
in the number of normalizable Majorana modes at each end is 
2. Here, the chain consists only of second-neighbor hopping 
and pairing. Note that there are now two extremal points of 
the bandstructure: one at fc = 0, the other at fc = tt. The 
winding number jumps from to 2 across this transition, and 
for A2 > 1 and Ai — 0, we see 2 normalizable Majorana zero 
energy modes at each end of the chain. 
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FIG. 5. (Color online)Topological phase transition across the point 
Ai — 1,X2 = 0. (a) Ai — 1^, A2 — 0. The quantities plotted are 
Efc (solid black line), /i (dashed black line), A(fc) (dashed blue line), 
and quasiparticle energy (solid red line) as a function of momentum 
in the one dimensional Brillouin zone, (b) The same quantities are 
plotted for Ai = 1^, A2 = 0. The associated Anderson pseudospin 
vector d{k) is drawn schematically below each plot. It is clear that 
in (a) the pseudospin does not wind along the Id Brillouin zone, i.e. 
W = 0, whereas in (b) it winds once, i.e. W = 1. 
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FIG. 6. (Color online)Topological phase transition in the vicinity of 
the point Ai = 0, A2 = 1. (a) A2 = 1", Ai = 0; (b) A2 = 1+, Ai = 
0. It is clear that in (a), the pseudospin does not wind along the 
Id Brillouin zone, i.e. = 0, whereas in (b) it winds twice, i.e. 
W = 2. 



Interestingly, our model hosts both BEC-BCS transitions 
and BEC-BCS crossovers. Only the former are topological 
transitions: these require that (1) an extremum of the band 
crosses the chemical potential, and (2) the pairing potential 
vanishes at the same momentum. If an extremum of the band 
crosses the chemical potential at a point where the gap does 
not vanish, the winding number will not change, since the total 
energy gap does not vanish. This is an example of a BCS-BEC 
crossover, and not a transition. This type of crossover is seen 
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FIG. 7. (Color online)Topological phase transition in the vicinity 
of the point Ai = 1.5, A2 = -1. (a) A2 = -l+,Ai = 1.5; (b) 
A2 — Ai = 1.5. In (a), the pseudospin does not wind along 

the Id Brillouin zone, i.e. = 0, whereas in (b) it winds twice, i.e. 
W = 2. 



near the line A2 = —1. An illustrative example is presented 
in Fig. [7] Here, the critical point Ai = 1.5 and A2 = —1 is 
studied. In Fig. |7^, the properties of the system are shown at 
Ai — 1.5 and A2 = —1^, where no normalizable Majorana 
zero mode occurs at the boundary. From the fact that Fermi 
points occur in this system, it is clear that the system is in the 
BCS regime. However, it is apparent from the form of the 
Anderson pseudospin that the winding number is identically 
zero. Thus, while this state is a BCS state, it is topologically 
equivalent to a BEC state which also has zero winding num- 
ber Thus, a crossover can connect this state to a BEC state. 
However, when Ai = 1.5 and A2 = — 1~, i.e. just below the 
critical point, we know from the analysis of previous sections 
that there are 2 normalizable Majorana fermion zero modes at 
each edge of the chain. This is also consistent with the wind- 
ing number of the Anderson pseudospin, which is T4^ = 2 
in this regime. We stress therefore, that a topological phase 
transition between 2 BCS states can occur. However so long 
as the bandstructure possesses inversion symmetry, it follows 
that such topological BCS-BCS transitions can only change 
the topological invariant by ±2. In a similar way, the critical 
line in Fig. ([TJ from b to d represents a topological transition 
across which the winding number changes by 2. Along this 
line, the gap vanishes at an incommensurate set of points in 
momentum space. In (c), the point Ai = 1 and A2 = —1 is 
considered. Here again, there are two band minima. How- 
ever, they occur at incommensurate momenta, ±kQ, symmet- 
ric about the origin. Therefore this critical point marks a tran- 
sition from to 2 Majorana zero modes at each end of the 
chain. This is the point (c) in Fig. dTli. As we approach the 
special multicritical point d in Fig. Jill (Ai = 2, A2 = — 1), 
these two incommensurate points move towards the origin. 
They meet at fc = 0, which is now a local maximum of the 
bandstructure. The 2 momenta ifcg meet at A; = at the 
multicritical point (Ai = 2, A2 = —1). In this way, all the 
topological phase transitions that occur in the model can be 



8 



(a) 



(0)8 

- i 


-2 





(b) I 
6 

4 

" 2 

-2 




1 


1 -1 

k/n 




k/it 


1 




(d) 8 
6 

m 4 
2 

-2 




1 


1 -1 

k/n 




k/n 


1 



FIG. 8. (Color online)Topological phase transition along the line 
Aa = -1. (a) Ai = 1.6, Aa = -1. (b) Ai = 1.8, A2 = -1. (c) 
Ai = 2.0, A2 = -1. (d) Ai = 2.2, A2 = -1. The associated 
Anderson pseudospin vector is not shown since it is not defined at 
points where the quasiparticle energy gap vanishes. However, the 
system in (d) is fully gapped and has W — 1, which is consistent 
with the analysis of previous sections. 



understood. 

Lastly, we study the nature of the special multi-critical point 
which occurs at (Ai = 2,A2 = —1). As this multicritical point 
is approached, the two incommensurate momenta which oc- 
cur around on either sides of k = 0, where the gap vanishes, 
approach each other at /c = 0. At the multicritical point, the 
momenta meet at fc = and annihilate each other, as shown 
in Fig. [8] Note that in Fig. [8](a-c), the Anderson pseudospin 
is not defined at points where the gap vanishes. Therefore, the 
winding number itself is not defined (this is consistent with 
the fact that this line represents a topological phase transi- 
tion). However, for Ai > 2, A2 = —1, the system is gapped 
everywhere, has a winding number W = 1, and possess one 
normalizable Majorana fermion zero mode at each end of the 
chain. 

To conclude, in this section, we have described a comple- 
mentary way in which the phase transitions in this model can 
be understood. Specifically, we have introduced the topologi- 
cal invariant corresponding to the winding number of the An- 
derson pseudospin vector around the Brillouin zone. Each of 
the phase boundaries of the spin model studied in this pa- 
per corresponds to regions where the winding number is ill 
defined, and across each critical point, the winding number 
changes by an integer The number of normalizable Majorana 
fermion zero energy modes localized to each end of the chain 
at a point in the AiA2-plane is exactly equal to the winding 
number at that point. We have also emphasized the point that 
a crossover can occur between a BEC and BCS system pro- 
vided that both have the same winding number [W ~ 0), and 
we have also demonstrated that there can be topological phase 
transitions from one type of BCS state to another 



C. Broken time reversal symmetry 

When the phase difference between the nearest neighbor 
and the second neighbor pairing amplitudes is non-zero, the 
BdG Hamiltonian takes the form 



H{k) 



4 

-afc - iPk 



-ak + iPk 
-a 



(34) 



where — A2sin6'sin2fc and j3k — Aisin/c + 

A2 cos0sin2A;. The particle-hole symmetry of the BdG 
Hamiltonian is 



aiH{k)ai = -H{-ky 



(35) 



The fc = 0, ±7r are special because they map onto themselves. 
Then from the equations 



= 0) = (1 - Ai - A2)a3, 
H{k = 'K) = (1 + Ai -A2)a3, 



(36) 
(37) 



it follows that the topological invariant is nfe=o tt sSii(?fe)- 
When ^0 = (1 - Ai - A2) and = (1 + Ai - A2) have 
opposite signs, we get n = 1, otherwise n — Q. As to the 
physical significance of fc 0, ±7r, it is similar to the case of 
topological insulators^'' where these points are termed "time 
reversal invariant" points. In a superconductor, k and -k states 
are paired, so the fc = 0, ±7r points are again special because 
they map onto themselves. This elegant topolo gical argument 



due to R. Roj0 confirms the results in Sec. 
Refs. mandEa 



IIIB 



see also 



V. CONCLUSIONS 

In this paper we have studied an exactly solvable spin 
Hamiltonian that is TFIM with an added three spin interac- 
tion. While such a spin interaction may appear to be artificial 
to the reader, such a term is generated in real-space renormal- 
ization group treatments of TFIM.''' Therefore, it is a physi- 
cally plausible coupling in a more reaUstic Hamiltonian and 
corresponds to longer ranged pairing and hybridization inter- 
actions among the fermions related to the spins via a Jordan- 
Wigner transformation. By analyzing the fermionized version 
of this spin Hamiltonian we have identified the quantum phase 
transitions in this system and to classify them according to the 
number of Majorana zero modes localized at each end. This 
number in turn is related to the winding number of the Ander- 
son Pseudospin unit vector along the one dimensional Bril- 
louin zone, so long as time-reversal symmetry for the spinless 
fermions is preserved. We have noted that when there are an 
even number of Majorana fermion modes at each end, there 
can be a crossover from a regime where the Majorana fermion 
wave-function decays in an oscillatory fashion (with an expo- 
nential envelope) to a regime where these modes decay expo- 
nentially without oscillation. Interestingly, at the crossover, 
the entanglement entropy vanishes identically. We stress that 
this crossover does not occur when there are an odd number 
of Majorana fermions at each end. Whether the vanishing 
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of the entanglement entropy is a necessary condition for this 
crossover remains to be understood. The degree to which such 
a crossover remains generic, or is ascribed to the integrability 
of the spin chain is unclear. Finally, an interesting possibility 
is that such crossovers may occur in higher dimensions in spin 
triplet superconductors in the presence of vortices, and other 
topological superconductors involving non-centrosymmetric 
systems. We shall relegate these studies to future work.. 

For unitary time reversal invariance, the topological argu- 
ment involving Anderson's pseudospin vector leads to wind- 
ing number Z. One might wonder if higher windings beyond 
n = 0, 1, 2 are possible as well. In principle, it is. To check, 
we added an even longer ranged term — Aacj^Ci+a + 
X^clcl^^ + h.c. Now, in addition to n = 0, 1, 2, we also get 
winding number n = 3 in appropriate regimes of the parame- 
ter space from explicit calculations of the BdG equation. It is 
quite likely that these higher order windings are energetically 
punished. The situation is very similar to the Xy-model in 
two dimensions for which higher order vorticity is suppressed 
by the chemical potential. Clearly phases with n = Z Ma- 
jorana zero modes are allowed by longer ranged Hamiltoni- 
ans. We find this phenomena intriguing, which deserves fur- 
ther attention. However, once the protection due to unitary 
time reversal invariance is removed the topological invariant 
collapses to Z2, with at most one Majorana zero mode at each 
end of an open chain. 

We have previously emphasized that while the solution of 
the spin model is exact, the fermonized version is a mean-field 
description of a p-wave superconductor whose exact solution 
requires treatment of fluctuation effects. In a recent paper it 



has been shown, however, that including fluctuation effects 
do not change the basic picture in a one-dimensional model.ESl 
Whether such a conclusion holds in higher dimensions, where 
Majorana zero modes are nucleated in the vortex cores of a 
Px + ipy superconductors, remains to be seen. We leave this 
problem for future research. 

An interesting question is whether or not the topologi- 
cal phases described here are perturbatively stable against 
weak interactions. We believe that they are, because they 
are gapped. In principle, for stronger interactions, the uni- 
tary time-reversal symmetry that protects the n = 2 phase 
can break spontaneously and destabilize it. The effect of 
stronger interactions in a specific model has been considered 
in Ref. 30 
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Appendix A: Broken time reversal invariance 



When there is a relative phase e between the nearest neighbor and the next-nearest neighbor pairing amphtudes the Majorana 
Hamiltonian is 



N 



N-1 



U = ^ hiQi + Ai ^ ha^+i + -^^^[{^ + cos 0)hi^iai+i 

1=1 i=l 1=2 

-(1 - cos6')ai_i&i+i + sin6'(ai_iai+i - 



(AI) 



Thus, we cannot simply set A2 to be complex in Eq. (I81. The Hamiltonian in the Majorana basis — 
(ai,6i,a2,62,a3,fe3,---)is 

...\ 



H 



( ^ 


1 








A2 
, 2 


sin^ 


-^(1 - cos 
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¥(1 


+ cos 9) 


-^sin6l 





-Ai 





1 














-1 







Ai 












-Ai 







1 














-1 






v 



(A2) 



To find the zero mode eigenvectors, we can try a solution of the form: Y^) = {Ai, Bi, A2, B2, A3, B3, • • • )^. However, the 
recursion relations turn out to be too complex to solve analytically. Thus, we resorted to numerical diagonalization of the BdG 
Hamiltonian in the main text. 



I 

Appendix B: Majorana zero modes in tlie presence of impurity for some special cases), while the regime with one Majorana 
In general, we find that when H of Eq. (Hh results in two 



Majorana zero mode, Himp in Eq. ( 22 1 destroys them (except 
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zero mode remains intact. 

Consider the general definition of the Majorana zero mode 
r = ^(Aiai+i3i6i), which is determined by the commutator 



— [-ffo + ^^imp, r] 

^2i^{Ai - AiA,+i - X2Ai+2)bi 
-2iBiai ~ 2i{B2 - XiBi)a2 (Bl) 
— 2i^(i?i+2 — XiBi+i — X2Bi)ai+2 
-\-2iX[Aj^mO,j ~ AjUj^jyi), 

which requires, in addition to the original recursion formula 

A, - AiA,+i - X2A,+2 = 0, (B2) 
B, - AiB,_i - X2B,_2 - 0, (i > J, m), (B3) 

new boundary conditions for B/s: 



Bi = for i < j, 



B, = XA 



Bj+m — XiBj^ni^i — A2-Bj+m-2 — ^^^j- 



(B4) 
(B5) 
(B6) 
(B7) 



(note that the Ai recursion relation is not affected by Bi's). 
Because of this change in the boundary conditions, we can no 
longer set Bi = for all i. Rather, for i > j + m, the general 
solution foiAj and Bj are of the form 



A, =C+qX + C-ql, 

B, ^c'+{i/q+y+cuyq-y 



(B8) 



where 1 — Xiq± — X2q±. — 0. 



We can now see how the impurity term Eq. ( 22 1 may de 



stroy the Majorana zero modes. Eq. ( |B8| l implies that if we 
had two Majorana zero modes without the impurity, which 
only requires \q±\ < 1, we will not have any normalizable 
Majorana zero mode due to the divergence of Bi unless we 
have Bj+jn-i = Bj+m = 0, which can occur only un- 
der special situations. On the other hand, if we had a sin- 
gle Majorana zero mode without the impurity, which means 
\q^\ < 1 < \q-\, the Majorana zero mode survives if 
Bi+m/Si+m-i = We have checked this explicitly for 

the special cases of m = 1, 2, 3. 



Appendix C: Impurity induced tunneling between two 
Majorana zero modes 

To consider the condition for the stability of two Majorana 
zero mode, we first not that, in the limit where the bulk gap 
is large, a semi-infinite chain can be regarded as a two-state 
system. This is because the two Majorana zero mode would 
form a single zero energy state, giving us energy degeneracy 
between the case where this zero energy state is occupied and 
the case where this zero energy state is vacant. Due to the 
fermion number parity conservation, perturbation cannot give 
rise to any off-diagonal term between the two states; all we 
can obtain is the energy difference between the occupied and 
vacant zero energy state. 

Therefore, an impurity term can annihilate the two Majo- 
rana zero modes if the mode expansion of this impurity term 
gives rise to dependence on the occupancy of the zero energy 
state. We know that, in absence of any impurity, the Hamilto- 
nian in Eq. Q gives us two Majorana zero modes near i — 1 
can be written down as the linear combination of only a^'s: 



(CI) 



where n = 1,2 and q £ TZ. Then, the annihilation operator 
of the zero energy state can be written as 

fo = (Ti + iT2)/2 = ^(ci, + ic2^)a,/2. (C2) 



What follows from this is that when we do the mode expan- 
sion on Majorana fermions on each site, only a/s receive con- 
tribution from the zero energy state whereas all bi's do not; 

CXD 

ai =(Cio/o + C*o/d) + (C«m/m + C*^fli)^ 

ra—l 

oo 

h^Y.^i[,Jm + i'Ul)- (C3) 

m— 1 

Any additional fermionic bilinear terms to the Hamiltonian 
cannot affect the zero modes unless mode expansion of such 
terms contain /q/q. 

We make a further restriction that we demand the fermionic 
bilinear to be local. The criterion for locality here is that, if 
our fermion operators are from sites i, j, they should satisfy 
\i-j\^0{\). 

This leads to the conclusion that only iaittj can gap out the 
zero modes, while iaibj and ibibj do not. (Conversely, if we 
had the right end of the semi-infinite chain, it is ibibj that gaps 
out the zero modes.) We see from Eq.(|C3|l 



zajflj = -i(c.,oc*o - c*iQC-jQ)flfo + (gapped) (C4) 

but the mode expansions of iaibj and ibibj do not have the 
/o/o term. 
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